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ESTIMATES  OF  BLOKVDOIVN  OF  QUASi-STATiC  PRESSURES 
IN  VENTED  CHAMBERS 

I.  INtRpDUCTldN 

In  earlier  work  on  estimating  quasi-static  pressures  within  vented  structure,  a dimensional 
analysis  was  made  and  scaled  curves  generated  ba'^ed  on  empirical  fits  to  data  for  paik  quasi- 
static pre^ures’and, blowdown  times  (see  Fig.  B5  and  B6  in  Ref.  'l).  Data  fromcReferences  2 
and  3 were  fitted  to  an  empirical  combination  of  scaled  parameters  which  seemed,  to  fit  the^ 
data  best. 

To  describe  more  completely  the  internal  pressures  applied  to  vented  structures  during 
the  blowdown  of  these  relatively  longterm  pressures  (long  term  compared  to  initial  reflected 
shock  pressures),  it  .would’be  desirable  to  be  able  to  estimate  the  time  histories,  including  their 
character  for  structures\which  were  better  vented -than  those  reported  In  References  2 and  3. 
Reference  3 reports  some  data  for  impulse  in  the  quasi-static  phase,  as  well  as  peak  quasi-static 
pressure  and  blowdown  times.-  These  authors  (Keenan  andTancreio)  also  note  that  they  db- 
tained'no  measurable  quasi-static  pressure  rises  for  values  of  (/lvem/^0  of  0.678  or  greater. 

No  measurements  have  apparently  been  made  in  the  range 

0.0775  < <0.678 

This  report  includes  additional  studies  on  blowdown  of  pressures  resulting  from  explosions 
In  vented  structures,  with  particular  emphasis  on  prediction  of  the  details  of  these  pressures  for 
design  of  suppressive  structures.  The  work  was  performed  for  the  Edgewood  Arsenal  suppressive 
structures  program  under  Contract  DAAA15-75-C-0083. 

II.  THEORY 

An  Analysis  is  available  in  the  literature  which  can  be  applied  with  some  modification  to 
this  problem.  Owczarek,'*  on  pp.  254-255  - considers  the  problem  of  prediction  of  pressure 
with  time  in  a tank  initially  filled  with  gas  at  high  prwsure,  and  then  emptied  by  the  sudden 
opening  of  a convergent  nozzle  in  the  tank.  wall.  He'^assumes*  a quasi-steady  flow  process,  a 
perfect  gas,  no  heat  transfer,  and  isentropic  expansion  in  the  hozzle.  Two  flow  regimes  apply, 
for  back  preMurep^.Iess  than,  and  equal/, )Or-greater  than,  critical  pressure  for  the  nozzle. 
Schefnatically„this  situation  is  shown  in  Figure  1 . 


*Hlsewhcre  in-Rcfercncc  '4,  Owczarck  aUo  considers  this  same  problem  with  shocks  and  rarefaction's  in  the  tank. 


i-IGURi;  I.  lOf-AliZED  VESSEL 

Owczarek's  dilutions  give  time  / as  a function  of  ratio  of,pn£ssu^at  this  time  to  outside 
ambient  pressure.f 


Po 
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These  solutions  are; 
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(4) 
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In  Eq.  (4),  is  the  time  at  which  P first  equals  the  upper  limit  in  Eq.  (^).  and  .v  is  a dimension- 
less group 

.V  = [/^>  " ‘ - l]  ''2  (6) 

tin  this  analysis,  all  pressures  aie  absolute  pressures. 
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In  these  equations.  /l,.-is  exit  area  of  the  nozzie.Tfhe  symbol  “a"  indicates  sound  velocity. 

7 is  ratio  of  specific  heats  for  gas  within  the  chamber,  is  the  universal  gas  constant,  and 
other  quantities  have  been  detined.  The  perfect  gas  law  requires  that 

</==  — = 7M  (7) 

P 

So  the  quantity 

(7/?.0|)'^"=«i  (8) 

i.e..  the  initial  sound  velocity  within'the  vessel.  The  set  of  Eqs.'(2)  and  (4)  applies  only  for 
initial  pressure  ratio  P|  ^eater  thaa  Usc- limit  in  Eq.  (5). 

With  some  as;,umptions.  adaptation  and  rcarrangement,  we  can  apply  Owczarek’s  cqua> 
tiohs  to  venting  of  suppa*ssive  structures.  Let  us  assume. that  the  nozzle  exit  area 

Ac  = pic  4s  (9) 

The  inlet  sound  velocity  is  an' important  parameter. in  Eqs.  (2)  and  (4).  In  suppressive 
structures,  this  is  no!  known  exactly  because  the  gas  being  vented  is  heated  air  mixed  with 
explosion  products.  But.  from  work  of  Proctor  and  Filleri^  it  is  apparent  that  the  great 
majority  of  the  ppjssure  rise  within  a vessefcontaining  ah  explosion  is  caused  by  addition  of 
the  heat  ener^  of  the  explosive  to  air,  rather  than  addition  of  ntols  of  explosive  products. 
Assuming  that, no  appreciable  venting  occurs  during  the  initial  stages  with  shock  reflections, 
this  process  represents  ^ constant  volume  heat  addition.  Such  a process  raises  the  absolute 
tempepture  of  a perfect  gas  in  direct  proportion  with  the  absolute  pre.ssure  rise.  From  Eq., 

( 3).  we  then  see  that 


Equations  (2)  and  (4)  can  be  rearpnged  in  dimensionless  form.  Let: 


no) 


T = 


fAcOi  _ /jaA  / Ac\ 

V " Vf'^V 


(ID 


where  f is  a dimensionless  time.  Then,  Eqs.  (2)  and  (4)  become 


2 /7+  IV^  +l)/|2<i  + i)| 

’^'"(7-1)  \~2~/ 


- I )f(2>) 

- 1 


fOTP>  1/ 


I) 


(2a) 
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and: 


^ ^ p<T  -n/(2t)  |q^4932_£(2v^  + 5K.v^  + I)'" 


'1^  | + '^*  + 0''*]} 


These  equations  express  dimensionless  time  as  a function  of  dimensionless  pressure  ijnd  ratio 
of  specific  hesits.  They  are  explicit  and  allow  solution  for  any  initial  pressure  ratio  P,  and  y. 
Plots  are  sK^owri  for  a range  of /’i  and  7 = 1.4  in  Figures  2 and  3.  From  these  figures,  one  can 
see  that  F = I for  nnife  values  of  r.  iie./biowdown  times  are  finite.  Scaled  durations  of  blow- 
down Tmax  be  obtaihed  from  Eq.  (2a)  and  (4a)  for  any  F|  by  determining'^  and  setting 
.V  = 0.  Values  for  scaled  duration  of  blowdown,  • and  time  to  reach  critical  pressure 
are  shown  in  Figure  4 as  funl'tions  of  F| . 


These.equations  change  the  definitiori-of  rin  Eq.  (11 ) to 


where  is  sound  velocity/in  the  outside  air,  /I,  is  internal  surface  area  of  the  structure,  and 
a,,  is  as  defined  in  Reference  i . From  Figure  4 and  our  previous  discussion,  we  note  that 
?ma\  is  a unique  function  of  F| . say 

?max=/i(F,;)  (1 3) 

However,  from  Eq.  (I2) 

^n.a,v=^  = /2(Fi)  (14) 


Furthennorc.  vre  can  define^q  simpler  Combination  of  parameters  as  a scaled  time.  i.e.. 

H^) 

so  that 


fnia\ 


(16) 
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!-IGURti4.  SCALED  BLOWDOWN  DURATIONS 


III.  COMPARISON  WITH  DATA 


The  function  Tm  ax  can  be  compared  with  measured  blowdown  times  for  vented  structures, 
as  from  Reference  3.  Theoritical  functional  forms  for  fmax,  are  plotted  versus  in  Figure  5, 
and  data  from  Reference  3 cast  into  this  form  are  also  plotted.:  Because  Referenced  cites  spread 
in  measurements,  we  can  indicate  this  spread  by  width  and  height  of  each  rectangle  in  Figure  S. 

The  experimental  data  describe  a different  functional  form  than  Owczarek’s  theory.  (Note  that 
the  vertical  scale  is  greatly  expanded,  and  that  deviation  from  the  theoretical  curve  is  always  with- 
in a factor  of  two.) 

Some  data  are  reported  in  Reference  S for  blowdown  times  for  cylindrical  vessels  with  end 
plates  having  various  sizes  of  holes.  These  investigators  ^ye  times  to  reach  critical  pressure,  rather 
than  atmospheric  pressure.  But,  we  can  easily  calculate  scaled  times  corresponding  to  this  pressure; 
they  are  r*  and  t„  defined  in  the  same  manner  as  in  Eqs.  ( 1 2)  and  ( 1 6)  respectively.  These  quan- 
tities! are  plotted  oti  Figures  4 and  5.  Scaled  data  for  t*  from  Reference  5 are  shown  on  Figure  5 
and  lje_generally  telpw  theiexperimental  valuesTonfmax. from^Reference  3. 

A much  more  approximate  analysis  of  blowdown  from  vented  structures  has  been  made  by 
Kinney  and  Sewell.^  They  assumed  that  there  was  no  regime  for  flow  transition,  and  venting: 
therefor^e  occurred  at  sonic  flow  velocity  down  to  outside  atmospheric  pressure.  They  also 
assumed  an  initial  source  temperature  independent  of  initial  pressure.  These  assumptions  yielded 
a simple  differential  equation  for  pressure  change, 


P^-C-P 


(17) 


Kinney  and  Sewell  integrated  this  equation  directly  and  obtained  the  solution'^ 

/n  P =fn  Pi  —CAt 

They  chose^v'  value  for  C = 2379  ft/sec,  giving  C = 2. 130.  Using  our  definition  for 

t = AT, 

Eq.  (18)  becomes 


(18) 


(19) 


^P  = APi~2A30r  (18a) 

and  it  yields  by  setting  P^  1.  This  value  is 

^niax  ~ 0*^695  A F|  (20) 

It  is  plotted  on  Figure  5 as  the  line  labeled  “Theory,  Kinney  and  Seweli”,  and  apparently  gives 
beeter  agreement  with  experiment  than  the  more  “exact”  theory  of  Owczarek. 

•We  have  converted  to  dimensionless  parameters,  w'lh  C = C/a0. 


r 


■<u. 


,'^cause  the  theories  of  Owczarek  or  Kinney  and  &well  predict  complete  time  histories 
of  blowdown  pr^ures,  these  histories  can  be  inte^ated  to  yield  predictions  of  gas  pressure 
impulses,  /g.  This  was  done  graphically  for  plots  siich  as  Figure  3 from  Owczarek’s  theory, 
and  analytically  for  the  simpler  equation  of  Kinney  and  Sewell.  The  latter  equation,  in 
dimensionless  form,  is 


7g=i  (e^>«-^  l)  -■ 


max 


(21) 


where 

(22, 

\ PoK  / 

Both  prediction  curves  are  shown  in  Figure  6,  together  with  points  from  Keenan  jind  Tancreto.^ 
Again,  the  very  approximate  formula  of  Kinney  and  Sewell  appears  to  fit  the  data  better  than 
the  more  “exact”  curve  based  on  Owczarek’s  equations. 


IV.  DISCUSSION 

From  the  analysis  presented  here,  it  seems  apparent  that  time  histories  of  pressures  within 
vented  explosion  containment  structures  can  be  scaled,  and  perhaps  can  be  estimated,  based  on 
a simple  gas  dynamic  analysis  by  Kinney  and  Sewell,*  or  a more  complex  analysis  by  Owozarek.** 
Tire  scaling  (dimensionless)  parameters  involved  are 


Owczarek’s  analysis  predicts  that 

?=/.  (24) 

and  gives  this  functional  dependence  (see  Figs.  2 and  3 for  a 7 = 1 .4).  Tire  combination  of 
scaled  parameters  in  the  second  term  in  the  parentheses  is  al.^o  a scaled  time. 
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FIGURE  6.  SCALED  GAS  PRESSURE  IMPULSE  VERSUS  SCALED  INITIAL  PRESSURE 


Because. P|  appears  iti  this  term,  a more  useful  scaled  combination  is  probably 


(12a) 


r 


= Al 


(14a) 


This  combination  also  appears  in  Kinneyaind  Sewell’s  simple  analysis. 

Both  prediction  equations  predict  finite  blowdown  durations,  w maxima  for  f and  f can 
be  easily  calculated.  These  maxima  and  scaled  times  to  reach  critical  pressure  are  shown  as 
functions  of  ?|  in  Figures  4 and  5 and  data  from  Reference  3 and  Reference  5 are  compared 
to  the  theoretical  predictions  in  Figure  S..  The  scaled  pressure-time  histories  can  be  integrated 
to  give  scaled  gas  impulse  Ig,  which  is  plotted  in  Figure  6 as  a function  of  Px  • This  parameter 
is  defined  as 


Limited  experimental  data  indicate  better  agreement  with  the  simple  equation  of  Kinney  and 
Sewell  than  with  the  more  exotic  solution  of  Owczarek.  Both,  of  course,  are  approximations 
to  the  complex  situation  occurring  during  venting. 


V.  CONCLUSIONS 

More  venting  pressure  data  should  soon  become  available  from  BRL  tests  in  the  Applied 
Technology  Program.  Also,  conversations  with  J.  Proctor  of  NSWC  Whjte  Oak  indicate  that 
there  is  a body  of  Norwegian  data  not  noted  here.  NSWC  is  also  looking  at  the  venting  process 
in  much  more  detail.  For  these  reasons,  curves  such  as  Figures  2 and  3 should  be  considered  as 
icntiiiive.  and  should  only  be  used  for  interim  estimates  of  venting  pressures,  gas  impulses,  and 
durations.  We  feel  that  the  primary  value  of  the  information  presented  in  this  rc'pprt  lies  in 
development  of  an  improved  method  of  presentation  of  venting  data  in  dimensionless  form, 
and  In  presentation  of  some  curves  and  formulas  which  give  preliminary  estimates  of  venting 
for  structures  with  large  vent  area  ratios. 
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